Quantum parallelism as a tool for ensemble spin dynamics calculations 
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Efficient simulations of quantum evolutions of spin- 1/2 systems are relevant for ensemble quantum 
computation as well as in typical NMR experiments. We propose an efficient method to calculate the 
dynamics of an observable provided that the initial excitation is "local" . It resorts a single entangled 
pure initial state built as a superposition, with random phases, of the pure elements that compose 
the mixture. This ensures self- averaging of any observable, drastically reducing the calculation time. 
The procedure is tested for two representative systems: a spin star (cluster with random long range 
interactions) and a spin ladder. 



PACS numbers: 03.67.Lx, 05.30.Ch, 75.40.Gb, 75.10.Pq 



One of the goals of quantum computers is the simula- 
tion of quantum systems. While quantum processing is 
ideally cast in terms of pure states, experimental realiza- 
tions are a major challenge [a met only fo r very small 
systems The alternative use of statistical mixtures 
of pure states, as the spin ensembles in standard NMR 
experiments d, 0, [^, 0], led to the development of the 
ensemble quantum computation (EQC) Q that allowed 
useful algorithm optimizations Recently, EQC was 

experimentally implemented in a 1 2-q ubits system [loj 
and much larger quantum registers were prepared to 
assess its stability against decoherence. Hence, an effi- 
cient evaluation of ensemble evolutions is needed. An- 
alytical solutions of ensemble dynamics, either from the 
integration of the Liouville von-Neumann equation [l^ 
or the alternative Keldysh formalism 1^ 14 1 are limited 
to special cases (e.g. [15|, ll6[). Thus, one has to resort 
to numerical solutions. One can describe an ensemble of 
M spins by using the 2*^ x 2*^ density matrix, but this 
quickly reaches a storage limit as M increases. Thus, a 
typical exact diagonalization method in a desktop com- 
puter slightly exceeds a dozen of qubits. Alternatively, 
the use of wave functions combined with the Trotter- 



Suzuki decomposition [17|, |18[ overcomes this limitation 
because it uses vectors of size 2^. However, the average 
over individual evolutions of a large number of compo- 
nents of the ensemble takes a long time. Here, this lim- 
itation is overcome by profiting of quantum parallelism 
[lit to evaluate the ensemble dynamics of any observable 
evolved from a "local" initial condition. The idea is that 
when evaluated on single pure states that are a superpo- 
sition of all the elements of the ensemble, these observ- 
ables become self-averaging. This reinforces the sugges- 
tions that one can avoid ensemble or thermal averages by 
using a single pure state 13, 21|. Ref. [13] considers a 
subsystem with the reduced density matrix derived from 
a pure state, where the subsystem is entangled with an 
environment which has a much bigger size. The resulting 
reduced density matrix describes the microcanonical en- 
semble without resorting to the equiprobability postulate 
of statistical mechanics. Following a similar inspiration. 



we focus on the non-equilibrium dynamics of any given 
observable. The key is that the initial non-equilibrium 
state has a "local" character, i.e. starting from an equi- 
librium ensemble state, a perturbing excitation acts on 
a small portion of the system. A number of useful tech- 
niques to simulate dynamics of a small system in the pres- 
ence of a mixed-bath that involve averages over random 
states [3 could be interpreted as particular cases. 
Specifically, an analytical justification of the numerically 
observed self-averaging property in some particular sys- 
tems [2^ follows from our results. 

While our method is general and can be applied to 
any mixed many-body system, we focus on the time evo- 
lution of the observable of greatest interest in NMR ex- 
periments, the local polarization ^, Ji, 5] . We consider 
two spin configurations that are representative of physical 
situations of contrasting topology. One is a spin ladder 
which is a variant of the linear chains which are exactly 
solvable 16 ,[3] and are known to have strong mesoscopic 
echoes The other is a spin star, a cluster with ran- 

dom long range interactions which is a fair representation 
of many molecular crystals [24|. Here, mesoscopic inter- 
ferences becomes less intense [25|. 

Ensemble vs. pure entangled state. — We take the en- 
semble of all the many-spin states l^*™) = |'0m) ® 
where m spins are in the state \ipm) and the are a 
base for the remaining M — m spins. These states have a 
statistical weight pt . The probability to find m' spins in 
the state \'4)m') at time t when the m spins were in state 
IV'm) at i = is 
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is the probability of finding 

to' spins in the state | ■(/;,„') within the many-spin state 

5*™ ^ at time t provided that, at < = 0, the to spins were 

in the state IV'm) within the state l^"™). The sum runs 
over all possible initial and final states. An example of 



2 



this is the local correlation function where m = m' = 1, 
\4'm) = |t)„ is the state of the n-th spin and I?/)™') — |t)„/ 
is the state of spin n! . The polarization of spin n! at time 
t provided that the spin n was up at time t = is given 
by P^n^i) = "^Wirii) - 1/2] il- The expression (P) 
involves D = 2^^^™ different dynamics for each of the 
initial states, see Fig. [Ija). This number is directly 
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FIG. 1; (Color online) Schemes of the quantum evolution of 
an ensemble (panel (a)) and a pure-state (panel (b)). Each 
l*r) = \i'm) ® contains a complete base, l^'i), of the 
M — m spins. 



related to the dimension of the Hilbert space. Thus, the 
number of computed evolutions increases exponentially 
with M. Our goal is to extract the same information 
in a shorter time. The parallelism implicit in quantum 
superpositions {T^ suggests that the desired correlation 
functions are contained in the dynamics of a single pure- 
state, see Fig. WO^)- The pure-state is built as an arbi- 
trary linear superposition of components of the ensemble. 



I.e. 
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-pure/ — J2iLi'^i\^T) where = ^p^e 
random ipi. Thus, the correlation function is given by 



with 
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Here {a} denotes the set of all the involved in the 
initial pure-state and D' — 2^"™ . Note that the sub- 
stantial difference between Eq. ([T]) and Eq. ([2]) is that 
the sum on i in the former is outside the square modulus 
while in the latter is inside. Rewriting Eq. ^ as 



f=ii=i 
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(3) 



the second term contains the initial correlations between 
the components of the initial state. These cross terms 
make the difference between W^')^(t) and wj^]^{t). By 
performing an extra average over Na realizations of the 
possible initial states, we obtain 



No. 



— yw^"^^ it). 
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Under this average, each element in the cross term 
goes to zero, hence W^"^(t) = hmAr^^ 



W^^?^ (t) 



Thus, by increasing the number Na of evolutions, ex- 
pression (HI) converges to W^,l^{t) as described by the 
central limit theorem. 
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The variance of wj^^l^{t) is 
Var = lim^„^^(W^r4W').. -(W^r^W)' • At this 



Na \ ™'™ -/A'„ 

stage, the four phases ipi summing up in each exponent 
are correlated. This enables terms where the exponent 
cancels out. These survive the average and contribute 



to Var 
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In the last expres- 
= Wi>i is the proba- 



where Pf ^ T,f ^/ /\*/ 

bility to find the state l^"™) after the subunitary evo- 
lution e^Hi/fi-p^Q-'^'Ht/h provided that the initial state 
was 1 5"™) . The projector Pj involves only part of the 
Hilbert space states. Hence, E?""^*'! — 1 implies 
Var = X]M'=iPjPi' Wei < EfPiniax{pj4 = max{pj/}. 
Thus, by using the Chebyshev's Inequality the prob- 



ability that 



VF^?^ (t) 



7V„ 



wens 
m'm 



(t) 



> e (i.e. the 



error exceeds a desired precision e) is lower than 
Var/ (iVaE^) < max{pi/} / (A^ae^) . For an homogeneous 
distribution Pi = l/2^^~™ hence max{pi/} / (iVaS^) = 
2/ {2^'^-"^Nae^) . The locality of the initial condition en- 
sures that M m, thus, as 2 
T^^?^ it) 



M-m 



TT/ens 
m'm V 



increases, one gets 
[t) in a single realization, i.e. the cross 
terms in Eq. ([3]) self-average to zero even for = 1. 

Spin systems with different coupling networks. — To il- 
lustrate the use of Eq. (|4]) we consider typical situations 
of high-field solid-state NMR. Here, the Hamiltonian is 
simplified by using a frame that eliminates the Zeeman 
contribution [l^]- We are left with the spin-spin interac- 



tion, H — J2 



M 
<3 
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, where 



bij/oij — represents an Ising-like coupling, Oij/bij = 
an XY Hamiltonian, Oij/bij = 1 the isotropic one, and 
Oij/bij — —2 a dipolar (secular) Hamiltonian truncated 
with respect to a Zeeman field along the z axis. The 
ensemble relevant for NMR experiments is in the infi- 
nite temperature limit [1, i-e., Pi = l/2^~'^ where 
{^i) are simple tensor product states in the Zeeman basis. 
The initial conditions are states with a local excitation 
at site n over a background level which is determined by 
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1, we have 



{a} 



M 

^prod/-iT>i0ni 

' 1=2 

A= (I j); + It); e^"^') with ipi random variables. Note that 



the zero magnetization of the other M —I spins [1, 0, Si ■ 

We calculate the local polarization of site n' at time t 
provided that it was polarized (n = n') at time t = in 
two different spin systems which have well differentiated 
kinds of dynamics: 

a) A ladder of spins interacting through an XY Hamil- 
tonian, as shown in Fig. [U^a). There, a^j = 0, hi^i+i — tion function (g]) is iWll°'^{t) 



)i , where |->>, 



- pure 



0i+M/2,i+M/2+l 



— 6t and b, 



i,i+M/2 



Here, the ex- 



act dynamics presents long lived recurrences (mesoscopic 
echoes) shown by the black line in Fig. [SJa), due to 
the high symmetry in the coupling topology [j, Q ■ The 
method also reproduces the exact solutions in isolated 
spin chains with both, pure XY or XY plus Ising interac- 
tions [i^l- Moreover, the results confirm that inclusion of 
Ising terms or interchain couplings leads to decoherence 
degrading the mesoscopic echoes 15| . 

b) A star system, see Fig. Hfb), in which all the 
spins interact with each other through a dipolar cou- 
pling Qij/bij = —2. The coupling intensities are given 
by a Gaussian random distribution with zero mean and 
variance a^. In this case, the local polarization decays 
with a rate proportional to the square root of the local 
second moment ctq = | (M — 1) 
nian and recurrences are negligible 
(exact solution) of Fig. [3l[b) shows the local polarization 
of this system. 



of the Hamilto- 
The black line 




FIG. 2: (Color online) Panel (a) shows the coupling network 
of a spin ladder. Panel (b) contains the coupling network of 
a spin star in which all the spins interact with each other. 



Testing the quantum parallelism. — In order to com- 
pare Eq. dH) to the ensemble average of Eq. ([1]), as 
well as its dependence on the choice of the phases ipi, 
we calculate the evolution for two types of initial states. 
Firstly, a pure entangled state is constructed by choos- 



ing ifi randomly. Thus, assuming n = 1, 



{"} \ 



pure , 



{a} 



-'v^ |t)^®|^r,) . The cor- 
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relation function, Eq. ([4]), calculated with this state 
is (WiT'(0>Ar„ giving the polarization (PiT(i)>Ar„ = 
2 ((W^ff (t))^^ - 1/2). The second case is a product (not 
entangled) state. It is built with the n-th spin up and all 
the others in a linear combination of spins up and down, 
with equal probability and arbitrary phase. Assuming 



this state can be rewritten in the form of 
the resulting phases (pi are correlated. Here, the correla- 

and the polarization 

- 1. 



is (Pf™'(i) 



= 2{Wir''{t) 



The local polarization, Pff'^(0 = 2 (VFff''(i) - 1/2) , 
obtained with Eq. ([!]), for the 14-spin ladder system is 
shown in Fig. with a black line. Square (red) and cir- 
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FIG. 3: (Color online) Local spin dynamics in a 14-spin sys- 
tem. The ensemble dynamics (solid line) is compared with 
that of entangled and product pure states. The square (red) 

and circle (green) scatter points correspond to {Pil°'^{t)^ 

and (Prr'(^))jY for = 1 respectively, (a) The extreme of 
a spin ladder with by/b^ ~ 1/10. The triangle scatter points 
(light magenta) correspond to (^Pii°'^{t)^ with A'^^ ~ 630, 

the lower value yielding the ensemble dynamics. A strong 
mesoscopic echo is evident, (b) A site in a spin star with 
random dipolar interactions. No mesoscopic echo is evident. 

cle (green) scatter lines correspond to the temporal evolu- 

jprod/ 



) be- tion of (^Pliit)J and (Pn*(*)>Ar„ respectively, with 
iV„ = 1. The agreement between Fff'(i) and (Pn'(^)>i 
is excellent, while (^Pii°'^ (t)'^ has a dynamics quite dif- 
ferent from that of the ensemble. The difference be- 
tween the dynamics of the two initial pure states is due 
to the different number of independent random phases 
of each state. In the random entangled state, there are 
2A/-1 independent phases that make the cancellation of 
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No, 



the second term in the rhs of Eq. ([3]) possible. How- 
ever, the number of independent phases for the prod- 
uct state is M — 1. This imphes that there arc multiple 
correlations between the phases in the cross terms in- 
hibiting their self cancellation. The triangle (light ma- 
genta) scatter line in Fig. El^a) shows the dynamics of 

It becomes indistinguishable from the ex- 
act dynamics provided that Na > 630. The relation be- 
tween NP/°'^ = 630 and iV™' = 1 is determined by the 
number of independent phases associated with the di- 
mension of the sampled portion of the Hilbert space [13] , 
i.e. 8192 = Ar™t2M-i iVP™d(M - 1) = 8190. 

Fig. mb) shows the local polarization for the spin 
star system. The complexity of this system washes out 
any possible recurrence for long times leading to a form 
of spin "diffusion". For 7V„ = 1 both (|pfr'^(t))^ and 

(Pff*(t))-^ are almost indistinguishable from the ensem- 
ble dynamics. This contrasts with the spin ladder where 

one would need (^Pfl°'^' 

Notably, (-Pi"*(0)i is an excelent approximant of the en- 
semble for both cases. This is because, in the star sys- 



PF^"^(t)) to get a fair description. 

630 



tem the cross terms 



of Eq. Q decay to a value of the order of 1/2^^ 
within a time scale determined by the Hamiltonian 
second moment ct^ = ^ctq. Thus, even the few M 
independent phases are enough to cancel the cross 
terms. In contrast, in the ladder system, the terms 



-mt/h 



present strong corre- 



lations and thus, the role of the phases becomes more 
relevant. 



In summary, we developed a method to overcome the 
limitations of the numerical calculations of an ensem- 
ble spin dynamics for large number of spins. Instead 
of evolving every one of the 2^~'^ initial states, when 
2^ ^ 2™, we evolve a single random entangled state. 
The procedure exploits the quantum parallelism implicit 
in quantum superpositions [l9| to reproduce the ensemble 
dynamics of any observable. This result supports a novel 
view of the foundation of equilibrium statistical mechan- 
ics [i^lllj. Moreover, even the non-equilibrium statisti- 
cal theory of the density matrix describing an ensemble 
in the thermodynamic limit could now be based on single 
states. Here, we observe that even for systems as small 
as 14 spins, the equivalence between a randomly corre- 
lated pure-state and an ensemble state holds. This is a 
consequence of the exponential increase of the dimension 
of the Hilbert space with the system size. The power of 
the method is enhanced when combined with the Trotter- 
Suzuki decomposition. We showed that the contribution 
of the extra correlations of the initial pure-state to the 
dynamics becomes negligible by increasing 2*^"™, the ra- 
tio between the size of the system Hilbert space and that 
of the subsystem where the non-equilibrium initial con- 



dition is supported. The method developed here allows 
for very efficient dynamical calculations of common ex- 
perimental situations where large ensembles are involved. 
Conversely, it prescribes possible pure input states for a 
quantum simulator to yield ensemble evolutions. 
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